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ABSTRACT 

Let K he a, 3-stranded knot (or link), and let L denote the number of crossings in K. Let ei and £2 be two 
positive real numbers such that £2 < 1- 

In this paper, we create two algorithms for computing the value of the Jones polynomial Vk (t) at all points 
t = exp (i(p) of the unit circle in the complex plane such that \(p\ < 2tt/3. 

The first algorithm, called the classical 3-stranded braid (3-SB) algorithm, is a classical deterministic 
algorithm that has time complexity O [L). The second, called the quantum 3-SB algorithm, is a quantum 
algorithm that computes an estimate of Vk (exp {i^p)) within a precision of £1 with a probability of success 
bounded below by 1 — £2- The execution time complexity of this algorithm is O (nL), where n is the ceiling 
function of (ln(4/£2)) /2£f. The compilation time complexity, i.e., an asymptotic measure of the amount 
of time to assemble the hardware that executes the algorithm, is O (L). 

1. INTRODUCTION 

Let if be a 3-stranded knot (or link), i.e., a knot formed by the closure 6 of a 3-stranded braid b, i.e., a braid 
b G B3. Let L be the length of the braid word b, i.e., the number of crossings in the knot (or link) K. Let £1 
and £2 be two positive real numbers such that £2 < 1. 

In this paper, we create two algorithms for computing the value of the Jones polynomial Vk {t) at all points 
t = e^"^ of the unit circle in the complex plane such that \lp\ < 

The first algorithm, called the classical 3-stranded braid (3-SB) algorithm, is a classical deterministic 
algorithm that has time complexity O (L). The second, called the quantum 3-SB algorithm, is a quantum 
algorithm that computes an estimate of Vk {'^^'^) within a precision of £1 with a probability of success bounded 
below by 1 — £2. The execution time complexity of this algorithm is O (nL), where n is the ceiling function 
of ^^^4^^^. The compilation time complexity, i.e., an asymptotic measure of the amount of time to assemble 
the hardware that executes the algorithm, is O (L). 
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2. THE BRAID GROUP 



The the n-stranded braid group B„ is the group generated by the symbols 

61, 62, • • • , bn-l 

subject to the following complete set of defining relations 

{bibj = bjbi for |i- j| > 1 

bibi+ibi = bi+ibibi+i for 1 < i < n 



This group can be described more informally in terms of diagrammatics as follows: Wc think of each braid 
as a hatbox with n black dots on top and another n on the bottom, and with each top black dot connected by 
a red string (i.e., a strand) to a bottom black dot. The strands are neither permitted to intersect nor to touch. 
Two such hatboxes (i.e., braids) are said to be equal if it is possible to continuously transform the strands of 
one braid into those of the other, without leaving the hatbox, without cutting and reconnecting the strands, and 
without permitting one strand to pass through or touch another. The product of two braids b and b' is defined 
by simply stacking the hatbox b on top of the hatbox b', thereby producing a new braid b ■ b'. Please refer to 
Figure 1. The generators fe,; are illustrated in Figure 2. Moreover, the defining relations for the braid group 
Bn are shown in Figures 3. The reader should take care to note that the hatbox is frequently not drawn, but is 
nonetheless understood to be there. 
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Figure 1. The product of two braids 



lil Hi - 111 

Figure 2. The generators of the n-stranded braid group B„ 
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Figure 3. A complete set of defining relations for the braid group -B„. 



Every braid b in the braid group i3„ can be written as a product of braid generators bi, 62, • ■ • , &n-i and 
their inverses b^^, 6^^, ... , 6^1 1, i.e., every braid b can be written in the form 



°j(l)°j(2) °jiL) 



where e (z) = ±1. We call such a product a braid word. 

Remark. We will later see that each such braid word can be thought of as a computer program which is to be 

compiled into an executable program. This resulting compiled program will in turn be executed to produce an 
approximation of the value of the Jones polynomial Jk (t) at a chosen point e*"^ on the unit circle. 

We define 

Definition 2.1. The writhe of a braid b, written Writhe{b), is defined as the sum of the exponents of a braid 
word representing the braid. In other words, 

Writhe =E^« 



For readers interested in learning more about the braid group, we refer the reader to Emil Artin's original 
defining papers, ^ ^ ^ as well as to the many books on braids and knot theory, such as for example.^ ^ 



3. HOW KNOTS AND BRAIDS ARE RELATED 



As one might suspect, knots and braids are very closely related to one another. 



Every braid b can be be transformed into a knot K by forming the closed braid b as shown in Figure 4. 




b b 

Figure 4. A braid b together with its closure b. 

This process can also be reversed. For Alexander developed a polytime algorithm for transforming an arbitrary 
knot K into a braid b having K as its closure. 

Theorem 3.1 (Alexander). Every knot (or link) is the closure of a braid. Such a braid can be found by 
a polynomial time algorithm 

Remark. Every gardener who neatly puts away his garden hose should no doubt be familiar with this algorithm. 

We should mention that it is possible that the closures of two different braids will produce the same knot. 
But this non-uniqueness is well understood. 

Theorem 3.2 (Markov). Two braids under braid closure produce the same knot (or link) if and only if one 
can be transformed into the other by applying a finite sequence of Markov moves 

We will not describe the Markov moves in this paper. For the reader interested in learning more about these 
moves, we suggest any one of the many books on knot theory.^^ 

4. THE TEMPERLEY-LIEB ALGEBRA 

Let d and A be indeterminate complex numbers such that d = —A"^ — A~^, and let 

Z[A,A-^] 

be the ring of Laurent polynomials with integer coefficients in the indeterminate A. Then the Temperley-Lieb 
algebra TLn (d) is the algebra with identity 1 over the Laurent ring Z [A, ^~^] generated by 

l,Ui,U2,...,Un-l 

subject to the following complete set of defining relations 

' UiUj = UjUi for \i-j\ > 1 

< UiUi±iUi = Ui 

^ Uf = dUi 

This algebra can be described more informally in much the same fashion as we did for the braid group: We 
think of the generators l,Ui,U2, . ■ . , Un-i as rectangles with n top and n bottom black dots, and with n disjoint 



red strings (i.e., strands) connecting distinct pairs of black points. The red strings are neither permitted to 
intersect nor to touch one another. However, they are now allowed to connect two top black dots or two bottom 
black dots, as well as connect a top black dot with a bottom black dot. The generators 1, Ui, U2, . . . , Un-i of 
the Temperley-Lieb algebra Tn{d) are shown in Figure 5. The reader should take care to note that the rectangle 
is frequently not drawn, but is nonetheless understood to be there. 
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Figure 5. The generators of the Temperley-Lieb algebra TL„(d). 



As wc did with braids, the product '•' of two such red stringed rectangles is defined simply by stacking 
one rectangle on top of another. However, unlike the braid group, there is one additional ingredient in the 
definition of the product. Each disjoint circle resulting from this process is removed from the rectangle, and 
replaced by multiplying the rectangle by the indeterminate d. In this way, we can construct all the red stringed 
boxes corresponding to all possible finite products of the generators 1, C/i, U2, ■ ■ ■ , Un-i- As before, two such red 
stringed rectangles are said to be equal if it is possible to continuously transform the strands of one rectangle 
into those of the other, without leaving the rectangle, without cutting and reconnecting the strands, and without 
letting one strand pass through another. Please refer to Figure 6. 
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Figure 6. Two examples of the product of Temperley-Lieb generators. 



Since TL„ {d) is an algebra, we also need to define what is meant by the sum '-|-' (linear combination) of two 
or more rectangles. This is done simply by formally writing down linear combinations of rectangles over the 
Laurent ring Z A"-'^], and then assuming that addition '+' distributes with respect to the product and 
that the scalar elements, i.e., the elements of the Laurent ring Z [A, , commute with all the rectangles and 
all the formal linear combinations of these rectangles. An example of one such linear combination is, 

{2A^ - 3A-^) l + + 7 A") Ui + {1 + A-'' - ^-1°) U1U2 , 

We should also mention that there exists a trace 

TrM:TLr,id)^Z[A,A-^] , 



called the Markov trace, from the Temperley-Lieb algebra TL„ (d) into the Laurent ring Z[A,A-'^]. This 
trace is defined by sending each rectangle to d*^"^, where k denotes the number of disjoint circles that occur 
when the closure of the rectangle is taken as indicated in Fig. 7. 



For readers interested in learning more about the Temperley-Lieb algebra TL„ (d), we refer them to the many 
books on knot theory, such as for example. 
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Figure 7. The Markov trace Tvm ■ TLn{d) — > Z [A, A''^] . 



5. THE JONES REPRESENTATION 

Vaughn Jones, using purely algebraic methods, constructed his Jones representation 

J:Bn^ Tin {d) 

of the braid group i3„ into the Temperley-Lieb algebra TL„ (n) by mapping each braid generator bi and its 
inverse into TL„ (d) as follows* 

hi ^ Al + A-Wi 

6^1 , — > A-^l + AU, 

He then used his representation J and the Markov trace Tvm to construct the Jones polynomial V {t) of a 
knot K (given by the closure 6 of a braid b) as 

V{t) = {-AY"''''^''^TrM{J{b)) 

where t = A~'^. 



'Actually to be perfectly correct, Jones wrote his original representation in a variable t which is related to our variable 
A by the equation t = A~^. 



Later, Kauffman created the now well known diagrammatic approach to the Temperley-Lieb algebra TL„ (d) 
and showed that his bracket polynomial (6) was intimately connected to the Jones polynomial via the formula 

{b)=TrM {J{b)) 



For readers interested in learning more about these topic, we refer them to the many books on knot theory, 
such as for example.^^ 



6. THE TEMPERLEY-LIEB ALGEBRA TL3 (D) 

We now describe a method for creating degree two representations of the Temperley-Lieb algebra TL3 (d) . These 
representation will in turn be used to create a unitary representation of the braid group -B3, and ultimately be 
used to construct a quantum algorithm for computing approximations of the values of the Jones polynomial on 
a large portion of the unit circle in the complex plane. 

From a previous section of this paper, we know that the 3 stranded Temperley-Lieb algebra TL3 (d) is 
generated by 

with the complete set of defining relations given by 

Uf=dUi and = dU2 

U1U2U1 = Ui and U2U1U2 = U2 
Moreover, the reader can verify the following proposition. 

Proposition 6.1. The elements 

l,UuU2,UiU2,U2Ul 

form a basis ofTLs {d) as a m,odule over the ring Z [A, . In other words, every element u) ofTLs (d) can 
be written as a linear combination of the form, 

CO = UIqI + COlUx + UJ2U2 + UJ\2UiU2 + L02lU2Ui = UqI + , 

where 

f-^O, Wi, lx'2, W12, U)21 

are uniquely determined elements of the ring Z [A, . 



7. A DEGREE 2 REPRESENTATION OF THE TEMPERLEY-LIEB ALGEBRA 

TL3 {D) 

We construct a degree 2 representation of the Temperley-Lieb algebra TL3 {d) as follows: 

Let |ei) and |e2) be non-orthogonal unit length vectors from a two dimensional Hilbert space Ti. From 
Schwartz's inequality, we immediately know that 



0<|(ei|e2)|<l 



Let S — i: |(ei|e2)| ^. It immediately follows that 

1 < 1^1 < oo 

Moreover, let a denote a complex number such that d = —o? — ot"^ . 



We temporarily digress to state a technical lemma that will be needed later in this paper. We leave the 
proof as an exercise for the reader. 

Lemma 7.1. Let S be a real number of magnitude \S\ > 1, and let a be a complex number such that S — —a^ — a~^. 
Then each of the following is a necessary and sufficient condition for a to lie on the unit circle: 

• 6 is a real number such that I < \6\ < 2. 

• There exist a 9 Cz [0, 27r] such that S ~ —2 cos (20). 

Thus, 

{a e C : 3 (5 such that \<\b\<1 andb = -o? - a"^} 
is equal to the following set of points on the unit circle 
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Also, as S ranges over all values such that 1 < |i5| < 2, a ^ ranges over two thirds of the unit circle, i.e., 
{ a-^ -.3 S such that 1 < \5\ < 2 and S = -a^ - a'^ } = | e'"^ : < y | 
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Figure 8. A plot of cos (261) for < 6* < 27r. 



We continue with the construction of our representation by using the unit length vectors |ei) and |e2) to 
create projection operators 

-El = |ei) (eil and £'2 = |e2) (62! 

These linear operators Ei and £'2 are elements of the endomorphism ring End {Ti.) = Mat (2, 2; C) of the Hilbert 
space Ti.. Since they are projection operators, they are Hermitian. By construction, they are of unit trace, i.e., 



tr{Ei) = l = tr{E2) 



where tr denotes the standard trace on End {H) = Mat (2, 2; C). The reader can also readily verify that 

tr{EiE2) = (5-2 = tr {E2E1) 

and that Ei and E2 satisfy the relations 

Ef = El and = E2 

E1E2E1 = 5 ^E\ and E2E1E2 = S ■^£'2 

It now follows that 

Theorem 7.2. Let 6 = ± |(ei|e2)|~"^ (hence, \6\ > 1), and let a be a complex number such that S = —o? — a~'^ . 
Then the map 

^a-TLsid) — > EndiH)^ Mat {2,2;C) 
Uj ^ SEj 

d I — > 6 
A I — > a 

is a well defined degree 2 representation of the Temperley-Lieb algebra TLs {d) . Moreover, we have 

tr($„(C/i)) = 5 = fr($„(;72)) 
and 

tr ($« {U1U2)) = l=tr ($« {U2U1)) 

Proposition 7.3. Let S = ±\{ei\e2)\~^ (hence, \S\ >1), and let a he a complex number such that 6 = —a^ —a~'^ . 
Moreover, let evala : Z \A, ^^^] — > C he the map defined by A 1 — > a. Then the diagram 

TL3 (d) ^ Mat (2, 2; C) 
TrM I Itr 
Z[A,A-^] C 

is almost commutative in the sense that, for each element lo e TL3 (d), 

evala oTru (w) =tr (w) + [5 —2) ljq, 

where denotes the coefficient of the generator 1 in w. 

8. A DEGREE 2 UNITARY REPRESENTATION OF THE THREE STRANDED 

BRAID GROUP B3 

In this section, we compose the above constructed representation with the Jones representation J to create a 
representation of the three stranded braid group B^. We then determine when this representation is unitary. 

We begin by quickly recalling that the 3-stranded braid group B3 is generated by the standard braid generators 

bi,b2 

with, in this case, the single defining relation 

616261 = 626162 



We also recall that the Jones representation 

Bs M TL3 id) 

is defined by 

bj I — > Al + A-^Uj 
bj^ ^ A-H + AUj 

where A is a indeterminate satisfying d = —A^ — A~'^. Thus, if we let 6 and a be as defined in the previous 
section, we have that 

'^-^ End{n) Mat {2,2;C) 
bj I — > al + a~^dEj 
bj I — > a~^I + adEj 

is a degree 2 representation of the braid group S3. Moreover, we have 



Proposition 8.1. Let J = ± |(ei|e2)| ^ (hence, \5\ > 1), and let a be a complex number such that S = —o^—a 
Then the degree 2 representation 

B3 "^-^ End (H) = Mat (2, 2; C) 

is a unitary representation of the braid group if and only if a lies on the unit circle in the complex plane. 
Proof. Since d is real and Ej is Hermitian, we have 

(al + a-^SEj)'' = al + a'^SEj 

So for unitarity, we must have 

al + oT^SEj = a"^7 + oi5Ej 
It now follows from the linear independence of /, Ei , E2 that $ o J is unitary if and only if 

a = 



□ 



From lemma 7.1, we have the following 



Corollary 8.2. Let 5 = ± |(ei|e2)| ^ (hence, \S\ > 1), and let a be a complex number such that 5 ■ 
Then the representation o J is unitary if and only if a = e*'', where 6 lies in the set 
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9. COMPUTING THE JONES POLYNOMIAL 



Recall that the Jones polynomial V (t) of a knot (or link) K given by the closure 6 of a braid word b is defined 



as 



where t = A~^. 

Thus, from Proposition 7.3, we know that the value of the Jones polynomial at a point t = e^"^ on the unit 
circle is given by 

y (eiV') = (_g3i.^^Ht/.e(6) ^^^j^^^ ^^^^ ^j^^^^ ^_^^,,^WritKeip) ^ ^ ^ ^ _ ^_^^^,^WritHeip) ^ 

where e'^ is a point on the unit circle such that e**' = (e*^) ^ = e~^'^. Prom lemma 1, we know that $e'« is 
only defined when |cos (20)| > i. Moreover, since y = —4^ mod 27r, it also follows from lemma 1 that ^^io is 
only defined when |<^| < 

Theorem 9.1. Let ip be a real number such that \(p\ < and let 6 be a real number such that (p = —AO mod 2-k. 
Let K be a knot (or link) given by the closure b of a 3-stranded braid b € B3. Then the value of the Jones 
polynomial V {t) for the knot (or link) K at t = e^"^ is given by 

V (e'^) = tr {U{b)) + {S-2) e'^^"'''<'>^ , 

where U is the unitary transformation 

U = U{b) = {^,.eoJ){b) . 

Let us now assume that \ip\ < ^ and that ip = —46 mod 2tt. Hence, U = U{b) is unitary. Thus, if the 
knot (or link) K is given by the closure 6 of a braid b defined by a word 
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k=l 



where 61, 62 are the generators of the braid group £3, and where e {k) = ±1 for fc = 1, 2, . . . , X, then the unitary 
transformation U = U (b) can be rewritten as 

fc=i 

where f/^-'^ denotes the unitary transformation (called an elementary gate) given by 

U^:>) = ($^,, o J) (bj) = ^^io {Al+A- lUj) = e'^I - 26"'^ cos (26') \ej) {ej\ 

= e'^I - 2e-'^ cos (261) Ej 



In summary, we have: 



Corollary 9.2. Let t = e^"^ be an arbitrary point on the unit circle in the complex plane. Let b be a 3-stranded 
braid (i.e., a braid b in B3) given by a braid word 

k=l 

and let K be the knot (or link) given by the closure b of the braid b. Then the value of the Jones polynomial 
V (t) of K att = e'*' is given by 

where U^^^ (j = 1,2) is the linear transformation 

U<^j) =e'^I- 2e-'^ cos (261) Ej , 

where I denotes the 2x2 identity matrix, and where Ej is the 2x2 Hermitian matrix \ej) {ej\. We have also 

shown that the linear transformations ?7^^\ ?7^^\ and U = 11^=1 (t/'"''''-'^)'^ are unitary if and only if \ip\ < 
When \ip\ < we will call U^^^ and U^"^^ elementary gates. 

Remark. Thus, the task of determining the value of the Jones polynomial at any point t = e^"^ such that < ^ 
reduces to the task of devising a quantum algorithm that computes the trace of the unitary transformation 

C/ = C/(6) = f[(?7(^W))'^'^ . 

k=l 



Corollary 9.3. Let K be a 3-stranded knot (or link), i.e., a knot (or link) given by the closure b of a 3- 
stranded braid b, i.e., a braid b G Bs. Then the formula found in the previous corollary gives a deterministic 
classical algorithm for computing the value of the Jones polynomial of K at all points of the unit circle 
in the complex plane of the form e*"^, where |<^| < Moreover, the time complexity of this algorithm is O (L), 
where L is the length of the word b, i.e., where L is the number of crossings in the knot (or link) K. We will 
call this algorithm the classical 3-stranded braid (3- SB) algorithm. 



10. TRACE ESTIMATION VIA THE HADAMARD TEST. 

In the past section, we hiWQ shown how to create the chissical 3-SB algorithm that computes the vahies of the 
Jones polynomial of a 3-stranded knot K on two thirds of the unit circle. In this section, we will show how to 
transform this classical algorithm into a corresponding quantum algorithm. 

We will now assume that \ip\ < ^ so that the elementary gates f/^^^and 1/^"^^ and also the gate U = U [b] = 

Y\k=i (U^^^''^^y^''^ are unitary. We know from the previous section that all we need to do to create a quantum 
3-SB algorithm is to devise a quantum procedure for estimating the trace trace (U) of the unitary transformation 
U. 



To accomplish this, we will use a trace estimation procedure called the Hadamard test. 



Let Ti be the two dimensional Hilbert space associated with the the unitary transformations U , and let 
{|/c) : fc = 0, 1} be a corresponding chosen orthonormal basis. Moreover, let JC denote the two dimensional Hilbert 
space associated with an ancillary qubit with chosen orthonormal basis {|0) , Then the trace estimation 

procedure, called the Hadamard test, is essentially defined by the two wiring diagrams found in Figures 9 and 
10. The wiring diagrams found in Figures 9 and 10 are two basic quantum algorithmic primitives for determining 
respectively the real part Re {trace{U)) and the imaginary part Im {trace{U)) of the trace trace [U] of U . The 
top qubit in each of these wiring diagrams denotes the ancillary qubit, and the bottom qubit \k) denotes a basis 
element of the Hilbert space Ti associated with U . (The top wire in each wiring diagram denotes the ancilla 
qubit.) Each box labeled by an 'i?' denotes the Hadamard gate 

The box labeled by an 'S" denotes the phase gate 

S = 

And finally, the controUed-t/ gate is given by the standard notation. 
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Figure 9. A quantum system for computing the real part of the diagonal element Ukk- 
Re (Ukk) = Prob [meas = 0) — Prob {meas = 1). 




Figure 10. A quantum system for computing the imaginary part of the diagonal element Ukk- 
Re {Ukk) — Prob {meas = 1) — Prob {meas = 0). 



The wiring diagram found in Figure 9 has been so designed as to compute the real part Re {Ukk) of the fc-th 
diagonal entry Ukk of U . For this wiring diagram has been so engineered that, when the output ancilla qubit is 
measured, then the resulting measured or 1 occurs with probability given by 

Prob{meas = Q) ^ i + i Re (fc|J7|fc) 

Prob{meas = l) = \-\Re{k\U\k) 

Thus, the difference of these two probabilities is the real part of the k-ih diagonal entry 



Prob {meas = 0) - Prob {meas = 1) = Re (fc|C/|fc) = Re {Ukk) 



If this procedure (i.e., preparation of the state |0) |A:), apphcation of the unitary transformation {H (g) l)-Contr — 
U- (if (g) 1), and measurement of the output ancilla qubit) is repeated n times, then the normaUzed number of 
O's minus the number of I's, i.e., 

#0's - #l's 
n 

becomes an ever better estimate of the real part Rc (Ukk) of the fc-th diagonal entry Ukk as the number of trials 
n becomes larger and larger. We will make this statement even more precise later. 



In like manner, the wiring diagram found in Figure 10 has been so designed to compute the imaginary part 
Im(Ukk) of the k-th diagonal entry Ukk of U. This wiring diagram has been engineered so that, if the output 
ancilla qubit is measured, then the resulting measured and 1 occur with probabilities given by 

Prob{meas = 0) = i - i Re (A:|i7|fc) 

Prob{meas=l) = ^ + ^Re{k\U\k) 

Thus, the difference of these two probabilities is the real part of the fc-th diagonal entry 

Prob {meas = 1) — Prob {meas = 0) = Im {k\U\k) = Im {Ukk) 

Much like as before, if this procedure (i.e., preparation of the state |0) application of the unitary transforma- 
tion (if (g) 1) • Contr — U ■ S- {H ^ 1), and measurement of the output ancilla qubit) is repeated n times, then 
the normalized number of I's minus the number of O's, i.e., 

#l's - #0's 

n 

becomes an ever better estimate of the imaginary part Im {Ukk) of the k-th diagonal entry Ukk as the number of 
trials n increases. 

We now focus entirely on the first wiring diagram, i.e.. Figure 9. But all that we will say can easily be 
rephrased for the second wiring diagram found in Figure 10. 

We continue by more formally reexpressing the wiring diagram of Figure 9 as the quantum subroutine 
QRE(7 ( k ) given below: 



Quantum Subroutine QRe;/ ( k ) 



Step 



Initialization 



IVo) = |0) \k) 



Step 1 Application of (g) / 



l^lj,) = {H <^ I) |^o) = _(|0)|fc) + |l)|A;)) 



Step 2 Application of Contr-J/ 



1^2) = (CONTR-[7) l^i) = (|0) \k) + |1) U \k)) 



Step 3 Application of H I 



1^3) = {H^I) 1^1) = |0) ( ^) \k) + |1) fi-^ ) \k) 



Step 4 Measure the ancilla qubit 



Resulting Measured Bit b 



Probability 



Resulting State \tp4 







i + iRe(Ukk) 



|0)(^)|k> 
^i + iRe{k|U|k> 



i-iRe(Ukk) 



|0)(^)|k) 
^i_iRe(k|Ulk> 



Step 5 Output the classical bit b and Stop 



Next we formalize the iteration procedure by defining the following quantum subroutine: 

Approx-Re-Trace;/ (n) 
LOOP k = 1..2 

Approx-Diag-Entry(/;;) = 
LOOP j = l..n 
b =QRe,7 (fc) 

Approx-Diag-Entry(/;;) = Approx-Diag-Entry(/;;) + (-l)* 

END loop j 
END loop k 

Output ( Approx-Diag-Entry(I) + Approx-Diag-Entry(2) ) 

End 

As mentioned earlier, quantum subroutines QlMu (k) and Approx-Im-Trace;/ (n) can be defined in a similar 
manor. 

We continue by recognizing that there is a certain amount of computational effort involved in creating the 
subroutine QREc/ { k ). For this, we need the following formal definition: 

Definition 10.1. The compilation time of a quantum algorithm is defined as the amount of time (computa- 
tional effort) required to assemble algorithm into hardware. 

Since the compilation time to assemble the gate U is asymptotically the number of elementary gates U^^^ in 
the product OLi (?7(^('=»)'^''\ we have 

Theorem 10.2. Let b be a 3-stranded braid, i.e., b G B3, and let K be the knot (or link) formed from the closure 
b of the braid b. Then the time complexity of compiling the braid word b into the quantum subroutine QRE[/ ( ) 
is O {L), where L is the length of the braid word b, i.e., where L is the number of crossings in the knot (or link) 



K. Moreover, the running time complexity of QRec/( ) is also 0{L). The same is true for the quantum 
subroutine QIm!7( ). 



Corollary 10.3. The quantum subroutine Approx-Re-TracE[/ (n) and Approx-Im-TracE[/ (n) are each 
of compile time complexity O {nL) and of run time complexity O [nL). 

Theorem 10.4. Let b be a 3-stranded braid, i.e., b G -B3, and let K be the knot (or link) formed from the closure 
b of the braid b. Let ei and £2 be to arbitrary chosen positive real numbers such that £2 < 1. Let n be an integer 
such that 

4 

Then with time complexity O (nL), the quantum algorithm Approx-Re-Tracej/ (n) will produce a random real 
number S^'^^ such that 

Probi^ 5^^^) -Re(trace(f/)) > ei ) < es 

In other words, the probability that Approx-Re-Trace;/ (n) will output a random real number Sn within ei 
of the real part Re {trace{U)) of the trace trace {U) is greater than 1 — €2- The same is true for the quantum 
subroutine Approx-Im-Traceu (n). 

Proof. Let Xi, X2, . . . , X„ be the n random variables corresponding to the n random output bits resulting 

from n independent executions of QRE[/ (1), and in like manner, let X^^i, X,i^2t ■ ■ ■ , be the n random 
variables corresponding to the n random output bits resulting from n independent executions of QRe^ (2). 

Thus, each of the first n random variables have the same probability pg^^ of being zero and the same probability 
p^^-* of being 1. In like manner, the last n of these random variables have the same probabilities p^'' and p^^-* of 
being or 1, respectively. Moreover, it is important to emphasize that the 2n random variables Xi, X2, . . . , X„, 
Xn+i,Xn+2, ■ ■ ■ 1 X2n are stochastically independent. 

The random variable associated with the random number 

#0's - #l's 



IS 



The reader can easily verify the mean fx of 5„ is given by 

M = Po^ - +Po^ - Pi^ = Re (C/11) + Re {U22) = Re {trace {U))) 
Prom Hoeffding's inequality, it follows that 

Prob ( \s^y - Re {trace{U))\ > ei ) < 2e-2(2^)'^?/(^f=i4) ^ ? 

Thus, when 



n > 



ln(2/e2) 



we have that 



Prob ( S^^""^ - Re {trace{U)) 



> ei < e2 



In like manner, a similar result can be proved for QlMf/. □ 



As a corollary, we have 



Corollary 10.5. Let b be a 3-stranded braid, i.e., b e B3, and let K be the knot (or link) formed from the 
closure b of the braid b. Let ei and £2 be to arbitrary chosen positive real numbers such that £2 < 1. Let n be 
an integer such that 

ln(4/£2) 



n > 



1e\ 



Then with time complexity O {nL), the quantum algorithms Approx-Re-TracE[/ (n) and Approx-Im-TracE[/ (n) 
will jointly produce random real numbers and Sn™^ such that 



Prob ( S^^''^ - Re {trace{U)) 



> £1 and 



S^""^ -lm{trace{U)) 



> €1 



< €2 



In other words, the probability that both Approx-Re-Trace;/ (n) and Approx-Im-TracE;/ (n) will output 
respectively random, real number sj^^'^^ and Sn'^^ within £1 of the real and imaginary parts of the trace trace (U) 
is greater than 1 — £2 • 



11. SUMMARY AND CONCLUSION 

Let be a 3-strandecl knot (or Hnk), i.e., a knot formed by the closure 6 of a 3-stranded braid b, i.e., a braid 
b & B^. Let L be the length of the braid word b, i.e., the number of crossings in the knot (or link) K. Let £1 
and £2 be two positive real numbers such that £2 < 1. 

Then in summary, wc have created two algorithms for computing the value of the Jones polynomial Vk (t) 
at all points t = e"^ of the unit circle in the complex plane such that < 

The first algorithm, called the classical 3-SB algorithm, is a classical deterministic algorithm that has time 
complexity O (L). The second, called the quantum 3-SB algorithm, is a quantum algorithm that computes 
an estimate of Vk (e"^) within a precision of £1 with a probability of success bounded below by 1 — £2. The 
execution time complexity of this algorithm is 0(nL), where n is the ceiling function of '"(^Z*^^) . The 
compilation time complexity, i.e., an asymptotic measure of the amount of time to assemble the hardware 
that executes the algorithm, is O {L). A pseudo code description of the quantum 3-stranded braid algorithm is 
given below. 

Quantum-3-SB-Algorithm(6, cp, e\, £2) 

Comment: b = braid word representing a 3-stranded braid s.t. K = b 

Comment: (p real number s.t. \ip\ < ^ 
Comment: £1 lower bound on the precision of the output 
Comment: £2 upper bound on the probability that the 
Comment: output is not within precision ci 

Comment: The output of this algorithm is with probability > 1 — £2 
Comment: a complex number within £1 of Vk (e"'') 

' ln(4/e2) 

U = Gate-Compile(6) 

Approx-Re-Trace;/ = Real-Part-Trace-Compile({7) 



ApPROX-lM-TRACEf/ = IMAGINARY-PART-TRACE-C0MPILE([/) 

AppROxReTr = Approx-Re-TracE{7 (n) 
AppROXImTr = Approx-Im-Tracec/ (n) 
W= Writhe {b) 
e = -if/A 
S = -2* cos (2 * e) 
ReExpS = cos (3 * 6* * W) 
ImExpS = sin (3 * 6* * W) 

ReJones =ReExp3*ApproxReTr-ImExp3 * ApproxImTr 
ReJones = (-1)^ * ^ ReJones + - 2) * cos {tp * W) 
ImJones =ImExp3*ApproxReTr+ReExp3 * ApproxImTr 
ImJones = (-1)^ * ^ ImJones - - 2) * sin {ip * W) 
Output( ReJones, ImJones ) 

End 
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